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Abstract 

We evaluate thermodynamical quantities and transport coefficients 
of a dense and hot hadronic matter based on an event generator 
URASiMA (Ultra-Relativistic AA colhsion Simulator based on Mul- 
tiple Scattering Algorithm). The statistical ensembles in equilibrium 
with fixed temperature and chemical potential are generated by impos- 
ing periodic boundary condition to the simulation of URASiMA, where 
energy density and baryon number density is conserved. Achievement 
of the thermal equilibrium and the chemical equilibrium are confirmed 
by the common value of slope parameter in the energy distributions 
and the saturation of the numbers of contained particles, respectively. 
By using the generated ensembles, we investigate the temperature 
dependence and the chemical potential dependence of the baryon dif- 
fusion constant of a dense and hot hadronic matter. 
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1 Introduction 



Physics of a high density and high temperature hadronic matter has been 
highly attracting in the context of both high energy nuclear collisions and 
cosmology as well as theoretical interest [|l| . In the recent ultra- relativistic 
nuclear collisions, though the main purpose should be confirmation of Quark- 
Gluon Plasma(QGP) state, physics of hot and/or dense hadronic state dom- 
inates the system. Hence, thermodynamical properties and transport coef- 
ficients of a hadronic matter are essentially important for the phenomeno- 
logical description of the space-time evolution of the produced exited region. 
In the cosmology, in addition to the global evolution of the early universe, 
baryon diffusion would play an important roll in the nucleosynthesis problem. 

Because of the highly non-perturbative property of a hot and dense hadronic 
state, investigation on the thermodynamical properties and transport coef- 
ficients has been hardly investigated. Numerical simulation based on Lat- 
tice gauge theory is a very powerful tool for the analysis of finite tem- 
perature QCD. Recently, transport coefficient of hot gluonic matter has 
been investigated 0. But even for the modern high-performance super- 
computer, lattice QCD evaluation of the transport coefficients of hadronic 
matter is very difficult, especially below T^. Furthermore, at finite density, 
present numerical scheme of lattice QCD is almost useless since inclusion 
of chemical potential makes lattice action complex although there are sev- 
eral new approaches have been proposed^j^]. In this paper, we evaluate 
the transport coefficients by using statistical ensembles generated by Ultra- 
Relativistic A-A collision simulator based on Multiple Scattering Algorithm 
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(URASiMA). Originally, URASiMA is an event generator for the nuclear col- 
lision experiments based on the Multi-Chain Model(MCM) of the hadrons[]5[. 
Some of us(N. S and O. M) has already discussed thermodynamical prop- 
erties of a hot-dense hadronic state based on a molecular dynamical sim- 
ulations of URASiMA with periodic condition[Q. Recently, some groups 
have been performed similar calculation with use of the different type of 
event-generator UrQMD0, where Hagedorn-type temperature saturation is 
reported. We improve URASiMA to recover detailed balance at temperature 
below two hundred MeV. As a result, Hagedorn-type behavior in the temper- 
ature disappears [^. This is the first calculation of the transport coefficient 
of a hot and dense hadronic matter based on an event generator. 

In section 2, we review URASiMA and explain how to make ensembles 
with finite density and finite temperature. Section 3 is devoted to the calcu- 
lation for nucleon diffusion constant through the first-kind fluctuation dissi- 
pation theorem. Section 4 is concluding remarks. 

2 URASiMA for Statistical Ensembles 

URASiMA is a relativistic event generator based on hadronic multi-chain 
model, which aims at describing nuclear-nuclear collision by the superpo- 
sition of hadronic collisions. Hadronic 2-body interactions are fundamen- 
tal building blocks of interactions in the model, and all parameters are so 
designed to reproduce experimental data of hadron-hadron collisions. Origi- 
nally, URASiMA contains 2-body process (2 incident particle and 2 out-going 
particles), decay process (1 incident particle and 2 out-going particles), reso- 
nance (2 incident particles and 1 out-going particle) and production process 
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(2 incident particles and n (> 3) out going particles). The production process 
is very important for the description of the multiple production at high ener- 
gies. On the other hand re-absorption processes ( n (> 3) incident particles 
and 2 out-going particles) thought to be unimportant in the collisions since 
system quickly expands and they have not been included in the simulation. 
On the other hand, in the generation of statistical ensembles in equilib- 
rium, detailed balance between processes is essentially important. Lack of 
re-absorption process leads one-way conversion of energy into particle pro- 
duction rather than heat-up. As a result, artificial temperature saturation 
occurs. 

Therefore, role of re-absorption processes is very important and we should 
take into account it. However exact inclusion of multi-particle re- absorption 
processes is very difficult. In order to treat them effectively, multi-particle 
productions and absorptions are treated as 2-body processes including res- 
onances with succeeding decays and/or preceding formations of the reso- 
nances. Here two body decay and formation of resonances are assumed. For 
example, A^A^ — > NNir is described as A^A^ NR followed by decay of 
R Ntt, where R denotes resonance. The reverse process of it is easily 
taken into account. In this approach, all the known inelastic cross-sections 
for baryon-baryon interactions up to ^/s < 3GeV, are reproduced. 

Table 1: Baryons, mesons and their resonances included in the URASiMA. 



nucleon 


A938 


A^1440 


A^1520 


A^1535 


A^1650 


A'l675 


A^1680 A'i720 


A 


^1232 


^1600 


■^1620 


^1700 


^1905 


^1910 


^1950 


meson 


TT 


V 




P770 
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For the higher energy, ^/s > 3GeV, in order to give appropriate total 
cross section, we need to take direct production process into account. Only 
this point, detailed balance is broken in our simulation, nevertheless, if tem- 
perature is much smaller than 3 GeV, the influence is negligibly small. For 
example , if the temperature of the system is 100 MeV , occurrence of such 
process is suppressed by factor of exp(— 30) and thus time scale to detect 
violation of detailed balance is very much longer than hadronic scale. 

In order to obtain equilibrium state, we put the system in a box and im- 
pose periodic condition to URASiMA as the space-like boundary condition. 
Initial distributions of particles are given by uniform random distribution of 
baryons in a phase space. Total energy and baryon number in the box are 
fixed at initial time and conserved through-out simulation. Though initial 
particles are only baryons, many masons are produced through interactions. 
After thermalization time-period about 100 fm/c, system seems to be station- 
ary. In order to confirm the achievement of equilibrium, we calculate energy 
distributions and particle numbers. Slope parameters of energy distribution 
of all particles become the same value in the accuracy of statistics(Fig. 1). 
Thus, we may call this value as the temperature of the system. The fact 
that numbers of species saturate indicates the achievement of chemical equi- 
librium(Fig. 2). Running URASiMA many times with the same total energy 
and total baryons in the box and taking the stationary configuration later 
than t = 150 fm/c, we obtain statistical ensemble with fixed temperature 
and fixed baryon number (chemical potential). 
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fig.l 



fig.2 



By using the ensembles obtained through above mentioned manner, we can 



3 Diffusion Constant 

According to the Kubo's Linear Response Theory, the correlation of the cur- 
rents stands for admittance of the system(first fluctuation dissipation the- 
orem) and equivalently, random-force correlation gives impedance(Second 
fluctuation dissipation theorem) 0. As the simplest example, we here focus 
our discussion to the diffusion constant. First fluctuation dissipation theorem 
tells us that diffusion constant D is given by current (velocity) correlation. 



evaluate thermodynamical quantities and equation of statesQ. 




(1) 



Average < • • • > is given by, 



1 
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(2) 



If the correlation decrease exponentially, i.e.. 



< v{t) ■ v{t + t') >oc exp ( ), 



(3) 
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with r being relaxation time, diffusion constant can be rewritten in the simple 
form, 

D = ^< vit) ■ vit) > T. (4) 
Usually, diffusion equation is given as, 

^J{t,x) = DV'f{t,x), (5) 

and diffusion constant D has dimension of [L'^/T]. Because of relativistic 
nature of our system, we should use /3 = 7- = ^ instead of v in eq.(|ip and 
D is obtained by, 



1 

D = - < I3{t) ■ I3{t + 1') > dt'c^. (6) 
3 JO 

= ^-<m-m>c'T. (7) 

^r(iii)'(ii)>- 

with c being the velocity of light. Figure 3 shows correlation function of the 
velocity of baryons. The figure indicates that exponential damping is very 
good approximation. Figure 4 displays the our results of baryon diffusion 
constant in a hot and dense hadronic matter. 



fig.3 



fig.4 
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Our results show clearer dependence on the baryon number density while 
dependence on energy density is mild. This result means importance of 
baryon-baryon collision process for the random walk of the baryons and thus 
non-linear diffusion process of baryons occurs. In this sense, we can state 
that baryon number density in our system is still high. In the inhomogeneous 
big-bang nucleosynthesis scenario, baryon-diffusion plays an important roll. 
The leading part of the scenario is played by the difference between proton 
diffusion and neutron diffusion ||T0|. In our simulation, strong interaction 
dominates the system and we assume charge independence in the strong 
interaction, hence, we can not discuss difference between proton and neutron. 
However obtained diffusion constant of baryon in our simulation can give 
some kind of restriction to the diffusion constants of both proton and neutron. 

From diffusion constant, we can calculate charge conduct ivity[jlT|. Figure 
5 shows baryon number conductivity ctb, 

= (9) 
KbT 

where tt-b is baryon number density, T is temperature and ks is Boltzmann 
constant (put as unity through out this paper), respectively. 



fig.5 



Therefore, if we want, we can discuss Joule heat and entropy production in 
the Baryonic circuit based on the above baryon number conductivity. 

Because fundamental system in URASiMA is high energy hadronic colli- 
sions, we use relativistic notations usually. However, diffusion equation (||) 
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is not Lorentz covariant and is available only on the special system i.e. local 
rest frame of the thermal medium. For the fuU-relativistic description of the 
space-time evolution of a hot and dense matter, we need to establish relativis- 
tic Navier-Stokes equation[^]. Taking correlation of appropriate currents, 
we can easily evaluate viscosities and heat conductivity in the same manner 



13 



4 Concluding Remarks 

Making use of statistical ensembles obtained by an event generator URASiMA, 
we evaluate diffusion constants of baryons in the hot and dense hadronic mat- 
ter. Our results show strong dependence on baryon number density and weak 
dependence on temperature. The temperature in our simulation is limited 
only small range, i.e., from 100 MeV to 200 MeV, and this fact can be one of 
the reasons why the change of diffusion constant of temperature is not clear. 
Strong baryon number density dependence indicates that, for the baryon 
diffusion process, baryon plays more important roll than light mesons. In 
this sense our simulation corresponds to high density region and non-linear 
diffusion process occurs. Calculation of the diffusion constants is the sim- 
plest examples of first fluctuation dissipation theorem. In principle, taking 
correlation of appropriate currents, i.e. energy flow, baryon number cur- 
rent, stress-tensor, etc., we can evaluate any kinds of transport coefficients. 
However, in relativistic transport theory, there exist several delicate points, 
e.g., relativistic property makes difference of mass and energy meaningless 
and, as a result, meaning of the "flow" of the fluid and "heat flow" become 
ambiguous |T2| [O]. The choice of the current depends on the macroscopic 
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phenomenological equations which contain transport coefficients. Once we 
estabhsh phenomenological equations for the high temperature and high den- 
sity hadronic matter, we can evaluate the appropriate transport coefficients 
in the same manner. Detailed discussion will be reported in our forthcoming 
paper I 
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Figure 1: The time evolution of the inverse slopes /3 ^ of Ng^s, A1232, P770 
and IT at ub = 0.1563 fm and etot = 0.938 GeV/fm . /S'^ is obtained 
by the energy distributions, ^ = 4^^^^^ = C exp(—PE). The dotted line 
stands for the /?~^ of vr. 
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Figure 2: The time evolution of number densities: (a) Etot 
0.313 [GeV/W] and (b) etot = 1-25 [GeV/fm% 
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Figure 3: Velocity correlation of the baryons as a function of time. Lines 
correspond to the fitted results by exponential function. Normalizations of 
the data are arbitrary. 
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Figure 4: Diffusion constant of baryons. 
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Figure 5: Baryon charge conductivity. 
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